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Net baryon number probability distribution near the chiral phase transition 
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We discuss the properties of the net baryon number probability distribution near the chiral phase 
transition. Our studies are performed within Landau theory, where the coefficients of the polynomial 
potential are parameterized so as to reproduce the mean-field and 0(4) scaling behaviors of the 
cumulants of the net baryon number. Within this model, we obtain the probability distribution 
for the net baryon number in a finite volume and study the effect of critical fluctuations on their 
properties. We show that in a mean-field approach, the singular part of the probability distribution 
can be obtained exactly and that it tends to broaden the distribution. By contrast, in the model 
with 0(4) scaling, the contribution of the singular part results in a narrowing of the net baryon 
number distribution. As expected, we find that the critical structures observed in the cumulants on 
the 0(4) critical line reflect modifications mainly in the tails of the distribution. We explore various 
methods for computing the probability distributions. 

PACS numbers: 25.75.Nq, 24.60.-k, 05.70.Jk 
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I. INTRODUCTION 

Fluctuations of conserved charges reflect the critical 
properties of a system. In particular, in a strongly in- 
teracting medium, fluctuations of the net baryon number 
and of the electric charge are valuable probes of the QCD 
phase transition [H-@|- Such fluctuations may provide a 
signature for the conjectured chiral critical point 043 , 
as well as for the residual criticality of the underlying 
0(4) transition dSEMl, expected in QCD at small 
densities in the limit of massless u and d quarks [ljl . 

Since fluctuations of conserved charges are experimen- 
tally accessible, they have been suggested as probes of 
the proximity of the chiral crossover to the freeze-out 
line in heavy- ion collisions [1(3, G3] ■ A particular role 
has been attributed to higher order cumulants, which ex- 
hibit an enhanced signature for criticality with increasing 
order 0, [HI [H|, [H| ■ In particular, the sixth- and higher- 
order cumulants of the net baryon number and of the 
electric charge vary rapidly in the transition region and 
even change sign close to the chiral crossover, already at 
vanishing baryon chemical potential 0, [IcJ O, H3] ■ 

Consequently, data on charge fluctuations in heavy- 
ion collisions at LHC and RHIC energies, may provide 
exp erimental evidence for the chiral transition of QCD [5, 
[l8[ . First measurements of fluctuations of the net baryon 
number, more precisely the net proton number, in heavy- 
ion collisions at RHIC have been obtained by the STAR 
Collaboration |19| . The data are, at least on a qualitative 
level, consistent with residual criticality, owing to the 
underlying chiral phase transition. 



Cumulants of conserved charges have also been studied 
theoretically. At small chemical potential fj,/T <C 1, they 
have been computed in first principle calculations, using 
lattice QCD [20l423| . while at large p., their prope rties 
were explored within effective chiral models [llT[l3l |2~I1 - 
[27j . which share the global symmetries of QCD. 

In statistical physics, the n-th order cumulant c n (T, /x) 
of a conserved charge N is obtained by differentiating 
the thermodynamic pressure p(T, fi) with respect to the 
corresponding chemical potential xx, 



d"[ P (?»/r 4 ] 
0(m/t)» 



(i) 



The pressure is related to the grand canonical partition 
function through p = (T/V)hiZ. The cumulants can 
also be expressed as polynomials in the central moments 
(((5 AO"), where SN = N— (N). The n-th moment (N n ) is 
linked to the net charge probability distribution through 



(N n ) = J2 N n P(N). 



(2) 



N=—oo 
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Cumulants of the net charge can exhibit singular prop- 
erties near a phase transition. Consequently, it is inter- 
esting to study the effect of critical fluctuations on the 
corresponding probability distributions. 

The aim of this paper is to explore the qualitative fea- 
tures of the probability distribution for the net baryon 
number near the chiral phase transition. Hence, we do 
not attempt the ambitious task to compute the probabil- 
ity distribution in QCD or in a QCD-likc effective model. 
Instead, we study the fluctuations within a more trans- 
parent framework, Landau theory of phase transitions. 
This approach does not provide a quantitative descrip- 
tion of the probability distribution. Nevertheless, the 
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qualitative aspects of our results are expected to be of order parameter is non-zero, a = ±-\/~ a(T, (J-) and the 
general validity. symmetry is spontaneously broken. The parametrization 

of a(T, /z) is chosen so that the broken symmetry phase 
is located below the critical line T = T c (fi), given by 
II. FLUCTUATIONS IN THE LANDAU a (T, /i) = 0. Consequently, the thermodynamic potential 

THEORY ; q 



In this work, we explore the influence of a second-order 
phase transition on the probability distribution of con- 
served charges within Landau theory of phase transitions. 

In Landau theory, the effective potential density uj at 
a given temperature T and baryon chemical potential fi 
is a polynomial in the order parameter cr, 



LUq = LU(T > T c (fi), (J,) = LU hg , 



(5) 



1 



w(T, n\ a) = cj bg + -a(T,/i)a + 



1 



(3) 



where uj^ g is a non-singular background contribution. Be- 
cause we focus on a regular second-order phase transition, 
far from a possible tri-critical point (TCP), we can trun- 
cate the polynomial in Eq. ([3]) after the fourth order term. 
In general, the effective quartic coupling can also depend 
on T and [i. However, near the critical line a(T c , /i c ) = 
and sufficiently far away from a TCP, such a dependence 
is irrelevant. Furthermore, at /x = we set the transi- 
tion temperature to T c = 0.15 GeV and parametrize the 
non-singular part of the Landau thermodynamic poten- 
tial density by 



co bg /T 4 = -2dcosh(/x/T). 



(4) 



For d — 7r 4 /30, the characteristic energy scale of the chi- 
ral phase transition in QCD is reproduced. This choice 
of parameters is clearly not unique. However, since we 
are interested in qualitative results, the precise choice of 
parameter values does not influence our conclusions. 
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FIG. 1: The phase diagram of the Landau model for small 
chemical potentials. The horizontal line indicates the path 
along which the cumulants and probability distributions of 
charge fluctuations are calculated. 

For positive a(T, fi), the minimum of the Landau po- 
tential ^ is located at a = 0, while for a(T, /z) < the 



uii = u(T < T c (p),n) = u hg - -|a(T,^)| 2 . (6) 

For the coefficient a(T,fi) we employ the following pa- 
rameterization 



o(T,/i) = Bign[d(I»] |d(T, M ) 



l-a/2 



where 



d(T,n) = - 



3-^-2cosh(/i/T) 



(7) 



(8) 



Note that for a = 0, a(T,fi) = d(T,fi). Moreover, for 
\T - T c \/T c <C 1 and [ijT c < 1, Eq. © reduces to 



tp, with A > and B > 0. 



d(T,fj.) ~ A{T-T c ) + Bn 2 
The scaling variable has been used frequently in the 
literature 0, 0] • 

The cosh(/i/T) terms in Eqs. {1} and dHJ account 
for the periodicity of the thermodynamic potential as 
a function of an imaginary baryon chemical potential, 
u)(T,i9T) with 9 = ih/T Q. We note that this period- 
icity is identical to that of the QCD thermodynamic po- 
tential 29;] , when \i is the baryon chemical potential. We 
do not consider the confinement transition, and identify 
/i with the chemical potential of the elementary fermion. 

The exponent a in Eq. ([7]) is introduced to parametrize 
the scaling properties of the thermodynamic potential. 
By tuning cv, we can switch between 0(4) and mean- 
field scaling 1 . For a = 0, the singular part of to exhibits 
mean-field scaling. On the other hand, for the 0(4) value 
of the critical exponent of the specific heat, a ~ —0.21, 
the Landau potential emulates the critical behavior of 
the 0(4) spin system in 3-dimensions, when the critical 
line is approached from the broken phase 2 . 

In Fig. [T]we show the critical line T c (/x c ) of the second 
order transition obtained in the Landau theory, using 
the parametrization ([7]) and © . The dependence of the 
first four cumulants on the chemical potential \i along a 
line of constant temperature is illustrated in Fig. [5] for 
T/T c = 0.98. 

Close to the critical point, at vanishing net baryon den- 
sity (|T - T c \/T c < 1 and /i/T c « 1), the singular part 



1 This admittedly ad hoc procedure provides a transparent frame- 
work for exploring the effect of critical scaling on the cumulants 
of the net baryon number. 

2 Since other critical exponents are not reproduced, this statement 
is limited to quantities governed by the critical exponent a. 
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FIG. 2: The baryon number density N/V and the first three cumulants c„ obtained from the Landau potential Eqs. ([5]) and 
© for T/T c — 0.98. The dashed lines shows the mean-field (MF) results with a — 0, while the full lines correspond to 0(4) 
scaling, with a = —0.21. 



of the thermodynamic potential (0 and (J6j) scales as 

1 



sing 



- i, 



|2-Q 



(9) 



for T < T c , while for T > T c the singular part vanishes. 
For /i ^ 0, the leading singularity of a given cumulant of 
the net baryon number is then given by 



(10) 



Using the 0(4) value for the exponent a ~ —0.21, the cu- 
mulants c„(T c , /i c ) diverge at the critical point for n > 3, 
as seen in Fig. [2J Moreover, in the model, the cumulants 
c n {T,n) are in general discontinuous at the transition 
point. 

With the mean- field critical exponent a = 0, the singu- 
lar parts of the cumulants remain finite near the critical 
point. Indeed, differentiating Eq. ^ for a — 0, one finds 
that for T < T c , the non- vanishing cumulants are given 
by 



sing _ 



•Bfi\t M \, d 



sing 



-B\t„ 



4 in s = QB 2 



2£?V, (11) 
(12) 



Consequently, also in a mean-field approach, the cumu- 
lants are in general discontinuous at the critical point, as 



seen in Fig. [5J but they remain finite on both sides of the 
transition. 

The singular behavior of the cumulants shown in Fig. [2] 
must be reflected in the corresponding probability distri- 
bution P(N) of the net baryon number. In the following, 
we compute the distributions P(N) in the mean-field and 
0(4) cases and explore the characteristic features, which 
are responsible for the critical behavior of the cumulants 
illustrated in Fig. [2] 



III. THE PROBABILITY DISTRIBUTION OF 
THE NET BARYON NUMBER 



Consider a sub- volume V of a grand canonical thermo- 
dynamic system of charged particles q and anti-particles 
q at a given temperature T and chemical potential fj,. The 
latter is related to the conserved net charge N = N q — Nq. 
The normalized probability distribution P(N) for finding 
a net charge N in the volume V is given in terms of the 
canonical Z(T, V,N) and grand canonical Z(T, V, (i) par- 



tition functions |30h32| 



P(N;T,n,V) 



Z(T, V, N)e^ N 
Z(T,V,ti) ' 



(13) 



where f3 = 1/T. 

The grand canonical and canonical partition functions 
are linked through the fugacity expansion, 



(14) 



N 



where A = e' 9 ^ is the fugacity parameter. 

As follows from Eq. (fT4|) . the canonical partition func- 
tion Z(T, V, N) is just the iV-th order coefficient in the 
Laurent expansion of Z(T, V, fi) about A = (/i = — oo). 
It follows that, given Z{T 1 V 1 /i), the canonical partition 
function Z(T, V, N) can be obtained using 



Z(T,V,N) = ^(f dX 



Z(T,V,n) 

X N+l ■ 



(15) 



where the integration contour C must lie in an annulus 
enclosing the origin in the complex A-plane. Inside the 
annulus the integrand Z(T, V, fi) must be analytic (see 
Fig. 13 

Taking C on the unit circle, A = e 10 , one finds the 
well-known result 130, 33], 



Z(T,V,N) 



1 

2^ 



271 



d9e 



-i8N 



Z{T,V,iT6). (16) 



We note that Eq. (fT6|) holds only when Z(T, V, fi) is ana- 
lytic on the unit circle. In general, the partition function 
in finite systems 3 , e.g. in lattice QCD simulations in a 
finite volume [35|-|4(|, is an analytic function of /i. Hence, 
in this case the integral (TToT) is well defined and indepen- 
dent of the radius of the contour. 

However, for the partition function Z = e~^ Voj , where 
lu is approximated by the thermodynamic potential den- 
sity in the thermodynamic limit, ([S])-®, while V is finite, 
there are singularities in the complex /i plane, which must 
be properly accounted for. This approximate treatment 
is justified for |iV| <C N* , where N* is a characteristic of 
the partition function of a finite system in the Yang-Lee 
theory of phase transitions 4 . 

To determine N* for a given system, a microscopic 
computation of the grand canonical partition function in 
a finite volume is needed. However, such a calculation 
cannot be carried out within the Landau model consid- 
ered here. Consequently, in the following we assume that 
N* is larger than the maximal N, needed in the numerical 
evaluation of the probability distribution P(N). We shall 
discuss the properties of the canonical partition function 
in this approximation. 



In a finite system, the partition function has Yang-Lee zeroes at 
complex im, which in the thermodynamic limit turn into cuts [44l . 
with branch points at the critical point and at \fi\ = oo. 
In the original work of C. N. Yang and T.D. Lee N* is the 
maximum number of classical particles that can be packed into 
the volume V. The Yang-Lee theory was later extended to a 
more general class of models. 



r 









FIG. 3: The integration contour C in Eq. (|15|) lies within a 
singularity free annulus in the complex A-plane, enclosing the 
origin. 



In the chiral limit, a critical point exists on the positive 
real fugacity axis |4a |. A = A c > 1. Charge conjugation 
symmetry of the partition function, /i — > —fx, implies 
that there is a corresponding critical point at A = 1/A C . 
This symmetry is respected also by the thermodynamic 
potential ©. 

At fixed temperature, the critical points at fj, = ±/x c , 
or equivalently at A = A^ 1 , are branch point singular- 
ities of the order parameter a in the complex /i (or A) 
plane (43[, with the critical exponent ft, a ~ (— . In 
the mean-field approach, this singularity is of the square 
root type, i.e. (3 — ^. Since in the mean- field approach 
^sing ~ W(T, [i)\ 2 9(—a(T, /j,)) and a ~ t^, the singular- 
ity of the thermodynamic potential at the critical point 
is a discontinuity in the higher derivatives, i.e. in the 
cumulants c n for n > 2, as shown in Fig. [21 

By contrast, for a = —0.21, the critical points are 
branch points of the thermodynamic potential, w S i ng ~ 
|£/j| 2 ~ Q - Consequently, the thermodynamic potential has 
cuts originating at the branch points. A convenient 
choice is to place the cuts on the real axis, between the 
critical point fj, = /i c and fi = oo as well as between 
fi = —fx c an d A* — — 00 • The corresponding cuts in the 
fugacity are located between A = A c > 1 and A = oo and 
between A = 1/A C and A = 0, respectively (42[. As a 
result, in the thermodynamic limit, the thermodynamic 
potentials of different phases (in the present case, the wo 
and tui) correspond to different Riemann surfaces con- 
nected through cuts [43l |. 

If the grand canonical partition function Z exhibits 
branch singularities, then the coefficients of the Laurent 
expansion Z(T, V, N), depend on the integration contour 
C in Eq. (j 15[) . In the broken-symmetry phase, the annu- 
lus is singularity free for pi = X c > 1 and p2 = 1/ A c < 1 
(cf. Fig. [3]). Thus, the coefficients of the Laurent expan- 
sion corresponding to the broken-symmetry phase, are 
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obtained using (IT6l) with Z — e~^ Vui1 . On the other 
hand, the partition function Z — e~^ v " of the symmet- 
ric phase is singularity free outside the annulus of the 
broken symmetry phase. Hence, the corresponding Lau- 
rent coefficients are obtained by integrating (|15[) along 
e.g. a circular contour with radius p < 1/A C or p > A. 

Thus, in general we obtain two competing partition 
functions for a given value of N . We note, however, that 
the partition function of a thermodynamic system must 
be unique. The reason for the ambiguity is the fact that 
we approximate u in Z by the grand canonical thermody- 
namic potential density in the thermodynamic limit. In 
the limit V — > oo, the uniqueness of the canonical parti- 
tion function Z is restored, since the one corresponding 
to a larger value of to is suppressed. 

In this paper, we consider only the canonical partition 
function obtained from the grand canonical one using 
(|15p for fi < p Cl i.e. for the broken symmetry phase. The 
relation between the probability distribution P{N) and 
the canonical partition function (|13p . implies that the 
structure of P(N) is entirely governed by the properties 
oiZ{T,V,N). 

The quantum statistics of fermion degrees of freedom, 
implies the existence of additional singularities in the 
complex [i plane, owing to the poles of the Fermi-Dirac 
distribution function [42| ■ The location of the "thermal" 
singularities may influence the possible choices of inte- 
gration contour in Eq. tjl 5|) , in particular for small values 
of the fermion mass. In the present model, which effec- 
tively corresponds to a Maxwell-Boltzmann treatment of 
the fermion degrees of freedom, these singularities are not 
accounted for. 



A. Probability distribution of the nonsingular 
potential 

The probability distribution corresponding to the non- 
singular part of the thermodynamic potential loq can 
be computed analytically. Indeed, using the generating 
function of the modified Bessel function I n {x), 



e 2 



n= — oo 



In(x)\ n , 



(17) 



one can directly expand the grand canonical partition 
function Z = e~^ Vu ° } without passing through the inte- 
gral representation (fl"5|) , and obtain [SO, l34j 



Z (T,V,N) = I N (2dVT 3 ). 



(18) 



Here we used JU) and cosh(/i/T) = (A + l/A)/2. 

The probability distribution from the nonsingular part 
of the Landau potential (|5|) is then obtained from Eq. (fT3)) 
as 

P NS (iV;T, V» = I N {2dVT 3 )e^ N+n °^ T , (19) 

where £! = Vljq — —2dVT 4 cosh(p/T) is the corre- 
sponding grand canonical potential. 



B. Probability distribution of the mean- field 
potential 

For the Landau potential ([6]) with the mean-field value 
for the exponent a = 0, the canonical partition function 
Z(T, V, N) may also be obtained analytically, using the 
same procedure as for the non-singular part. 

The partition function of the singular case is written 
in the following from 



exp <^ VT 



d A + 



-a 2 (7» 



(20) 



By expanding the argument of the exponential (here we 
put t = 1 - T/T c ) 



«(T.//)|- = (/ + 2)- -2- 2(1 2) j A + I J + A 2 + 1 



A 2 ' 
(21) 



and using Eqs. (|T7|) and (fT6)l . we find 



(T,V,N) = e^ t+2 y 2+ ^ 



J2 I N -2i[{2d-t-2)VT 3 }I ( 



{VI* 



V 2 



(22) 



The probability distribution P MF (N) can be then com- 
puted using Eqs. (|2T)|) . (122]) and (fT3)) . 



P MF (N;T,V,p) = e VT3 ( t -^- d )( x+ i)-H^+^) 2 +^ 

-y T 3 



I N -2i[(2d-t-2)VT 3 ]I e 



(23) 



In Fig. |4] we show the resulting probability distributions 
P MF (iV) and P NS (A0 a t M = for V = 30 fm 3 . 

The mean-field probability distribution P MF (N) is 
clearly broader than P NS (iV). This feature is also evi- 
dent in the ratio P MF (N) / 'P NS (AT) , shown in the right 
panel of Fig. |U However, for small N, \N\ < 30, the 
trend is reversed and P is narrower than P NS (7V). 



C. Probability distribution in the critical case 

We have shown, that for the regular as well as singular 
part of the potential with the mean-field exponent a = 0, 
the net charge probability distributions can be obtained 
in closed form. The singular part with a = —0.21, how- 
ever, does not allow for an analytic calculations. In this 
case, we can apply a numerical or some approximate ana- 
lytical method to obtain the canonical partition function 
and the corresponding probability distribution P a (N). 
The most direct approach to compute Z(T, V, N), is the 
numerical integration of Eq. (|16[) . There are, however, 
limitations of this method, in particular at large N, due 
to the oscillatory structure of the integrand. We adopt 



6 



1 

1e-10 

1e-20 

% 1e-30 
a. 

1e-40 
1e-50 
1e-60 



H=0 

V=30 fm 3 
MF 

Non-singular 



S 



-150 -100 -50 




N 



50 100 150 



1.3 
1.25 
£ 1-2 

g 1.1 

1L 

Q- 1.05 
1 

0.95 



, MF --- 
i Non-singular 


/ 

r 
( 


n=o 

. \ V=30fm 3 

\ 

X 

\ 
\ 

\ 

\ 

\ 

\ 

N ^ — 


; 

/ ■ 
/ 

; ■ 

/ 
/ 
/ 

/ 

/ 

/ 

/ 

/ 

/ 



-40 



-20 




N 



20 



40 



FIG. 4: Left : Probability distributions of the net baryon number in the Landau model. The left-hand figure shows the 
singular contribution in the mean-field case P MF (N) from Eq. fl23j (dashed-line) and the nonsingular P NS (iV) from Eq. |19]) 
(dotted-line). The right-hand figure shows the ratio of these probability distributions. The calculations were done at /j, = 
and for sub- volume V = 30 fm 3 . 
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FIG. 5: Comparison of the probability distributions of the net 
baryon number P(N) in the Landau model with the mean- 
field exponent calculated through different methods, see text. 
Each P(N) is normalized to an exact result P MF (N) fr om 
Eq. (ESI). 



FIG. 6: The probability distribution calculated by the nu- 
merical integration from Eq. (|16|) in Landau model with 
a = —0.21 and normalized to the non-singular probability 
distribution P NS (JV) from Eq. (fT9)l . 



the qawo subroutine in QUADPACK package library for 
computing oscillatory integrals. 

To test the numerical method, we first compare the 
probability distribution P MF (N) obtained analytically 
in Eq. (|23|) with the numerical integration of Eq. (fl~6|) . 
In Fig. [5] we show the ratio of the resulting probability 
distributions. The numerical integration reproduces the 
analytical result up to \N\ ~ 60, but fails for larger N. 
Thus, we conclude that the numerical integration can 
be used with confidence also to compute P a (N) for the 
model with a ~ —0.21 up to this value of N. 

In Fig. [6] we show the probability distribution P a (N) 
obtained through numerical integration of Eq. (1161) . 
normalized to the non-sing ular distribution P NS (A), 
Eq. (|2U)) . This ratio exhibits quite different behavior 



from that seen in Fig. 0] for the mean- field model. While 
P MF (N) keeps increasing compared to P NS (iV), the sin- 
gular distribution P a (N) decreases relative to the non- 
singular background for large N. This indicates, that for 
the critical exponent a ~ —0.21, the singular probability 
distributions is narrower than the non-singular reference 
distribution. 

The canonical partition function Z(T, V, N) can also 
be computed by applying the method of steepest descent 
to the integral in Eq. (|15j) . To this end, it is convenient 
to change the integration from the complex A- to the 
complex /i-plane. The integration along the closed con- 
tour C, characterized by p 2 < |A| < p\ (see Fig. [3] ), 
is transformed into a line integral from \i = /ir — inT 
to n = jiR + z7rT, where /ii? is restricted to the range 
— /j, c < fiR < fi c (see Fig. [?])■ The canonical partition 
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FIG. 7: The integration paths in the complex /i-plane of the integral (124[l . The left-hand figure, shows the path in the (MF) 
case where a = in (|24[) . The solid line along the vertical axis denotes the original integration line and the shifted solid curve 
is the steepest descent path for the saddle point y, s on the real /i axis. The right hand figure, shows the path for a = —0.21 
in (|24[) . The lines on the horizontal axis, starting from fi = ±/i c to ±oo, are the cuts of the integrand. The saddle points for 
| TV" | > N c are indicated by the crosses. The steepest descent path at the saddle point is indicated by the arrow. 
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FIG. 8: The real part (top) and imaginary part (bottom) of 
the saddle point ji s in the complex /i-plain as a function of N 
for the mean-field and the critical case with a = —0.21. 



function is then given by 

Z(T,V,N) = JLf c ^e 
where 

f{fi) = -(a; 1 (T )/ x)+ At JV/y)/T 4 



vt 3 /(m) 



(24) 



2dcoM»/T) + \\a\ 2 {T,ri-^L. (25) 



For sufficiently large values of the pre- factor VT 3 , the 
integration ([24]) can be evaluated using the saddle-point 
approximation. For V = 30 fm 3 and T = 150 MeV, one 
finds VT 3 ~ 13. This is reasonably large, so that we 
expect acceptable results in the saddle point approxima- 
tion. 

Let fig — fi s (T, V, N) be the solution of the saddle point 
condition 



dm 



dfx 



= 0. 



Then, the canonical partition function reads 

1 



Z(T,V,N) = 



y/2nVT^\f"(n a )\ 



(26) 



(27) 



Equation (f25j) . implies that the saddle point condition 
(|26| is given by 



N = - ^ 



(XT 

df 1 / T,V 



(28) 



the thermodynamic relation for the average baryon num- 
ber N — (N) in the grand canonical ensemble. For 
/i < /i c the relevant thermodynamic potential is that of 
the ordered phase, ft = Qi. Thus, for N < N c , where 
N c ~ 11 is dehned by the baryon number at the tran- 
sition point /i = He, the location of the saddle point /j, s 
is given by the density, shown in Fig. [2] On the other 
hand, for N > N c the saddle point condition is modified, 
as discussed below. 

In the mean-field parametrization we find that the sad- 
dle point on the ill Riemann sheet is always located on 
the real axis, reflecting the lack of branch point singu- 
larities of the thermodynamic potential (see Fig. [7]). On 
the other hand, in the critical case fi s is real for N < N c , 
but for N > N c there are two complex conjugate saddle 
points on either side of the cut, as shown in Fig. [7J The 
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imaginary part of the saddle points fi s as a function of 
N is shown in Fig. [5] in the mean- field and in the critical 
case. 

To test the validity, and the range of applicability of 
the saddle point method, we apply this approach to com- 
pute the canonical partition function and the correspond- 
ing probability distribution for the mean-field model. In 
Fig. [5] we show a comparison of P(N) calculated using 
the saddle point integration with the analytical result 
P MF (iV), Eq. ([23]). The accuracy of the saddle point in- 
tegration decreases with increasing |iV| . With the param- 
eters used in our calculations, the saddle point approxi- 
mation reproduces P MF (N) with sufficient accuracy, only 
for \N\ < 20. The range of applicability in N is much 
smaller than that of direct numerical integration. We 
have verified, that by increasing VT 3 , the range of appli- 
cability of the saddle-point approximation grows to larger 
\N\. 

We note that the probability distribution P a (N) ob- 
tained with the saddle point method exhibits larger devi- 
ations from the numerical results than in the mean-field 
case. The main problem is that in this case, for N > N c 
there are two complex conjugate saddle points, located 
close to the cut, as illustrated in Fig. [5J Consequently, 
there is also, in general a contribution to the integral 
emanating from the cut, which is neglected in the saddle 
point approximation. 

Neglecting the cut and summing up the contributions 
from the two saddle points, the resulting partition func- 
tion OH) reads, 



Z(T,V,N) 



v /2nVT 5 \f"(^ s )\ 
x cos(VT 3 Im[/(/i s )] 



,VT 3 Re[/(^)] 



(29) 



The cos(a:) term in the above expression yields an oscil- 
latory contribution to P(N). This implies, that in some 
range of parameters, P{N) is negative. This indicates, 
that the saddle point approximation in the present form 
has a rather restricted range of applicability. 

In general, the appearance of the oscillatory factor in 
P(N) is attributed to the contribution of the complex 
saddle points, resulting from the non-analyticity of the 
thermodynamic potential due to a phase transition. It 
is expected that in the physical case, this problem would 
be removed when the contribution of the cut is included. 

In the present case, however, the onset of negative 
probabilities cannot be attributed solely to the complex 
structure of the saddle points. In fact, in our model, 
the negative probability appears for \N\ > 50 also when 
the canonical partition function is evaluated numerically. 
Thus, the unphysical behavior of the canonical partition 
function seen in Eq. ([29)) might not only be attributed to 
the saddle point approximation but possibly also to an 
artefact of our model for the grand canonical partition 
function. 
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FIG. 9: Cumulants calculated from the probability distribu- 
tion P a (N) in Landau model for a = -0.21. The P a (N) was 
calculated numerically up to \N\ = 50. The curve labelled 
GC is obtained directly from the grand canonical partition 
function using Eq. {T}. 



D. Cumulants calculated from P(N) 

In the preceding section we have derived the probabil- 
ity distributions for the Landau potential applying two 
different methods, the numerical integration and the sad- 
dle point approximation. In general, the probability dis- 
tribution exhibits characteristic structures near the phase 
transition, which are reflected in the critical behaviour 
of the cumulants. To explore the role of the critical fluc- 
tuations, we compute the cumulants directly from the 
probability distributions P(N), 



c 4 (T, n) 



1 
1 
1 



((<5iV) 2 ), (30) 
((SNf), (31) 
yT3 [{(6N) 4 )-3((6N) 2 ) 2 ], (32) 



where the average ((SN) n ) is defined in Eq. ([2]). 

The analytic form of the mean-field probability distri- 
bution P MF (N), reproduces the properties of cumulants 
shown in Fig. [3] The approximate P MF (N) obtained 
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by numerical integration, yields cumulants c„ , which are 
consistent with exact results except for the immediate 
neighbourhood of the phase transition. Furthermore, the 
error of the saddle point approximation is at the few per- 
cent level for small /i, but increases as /i c is approached. 
Such differences are expected, since the numerical and 
saddle point probability distributions do not reproduce 
the analytic results at large N. 

In the model with the 0(4) critical exponent a = 
—0.21, deviations are more pronounced. This is illus- 
trated in Fig. [31 where the second- and the fourth-order 
cumulant, obtained using Eqs. (|3D ]) -([52" j) with \N\ < 50 
and the probability distribution evaluated numerically, 
are confronted with the exact results, computed using 
Eq. ((TJ with p = —oJi. The second order cumulant, 
which remains finite at the 0(4) critical point, is well 
reproduced, while the for fourth order cumulant, which 
diverges at the transition point at finite density, we find 
agreement only at small \x < /i c . Thus, the approximate 
cumulant completely fails to reproduce the critical struc- 
ture near /z c . This is closely related to the fact, that the 
P(N) obtained by numerical integration does not repro- 
duce the large N behaviour, as illustrated in Fig. [5] for 
the mean-field model. Since the fugacity factor e^ N l T in 
Eq. (|13|) enhances the large N part of P(N) for /x > 0, 
P(N) at \N\ > 50, which is negligible in P(N) at fx <C (j, C) 
becomes necessary to reproduce the higher order cumu- 
lant. 

Moreover, in our calculation the sub-volume is rather 
small and therefore the critical structure of the higher- 
order cumulants at the second order transition is not 
clearly exhibited. In order to reproduce the singular 
structure at the critical point in the cumulants of P(N), 
the probability distribution in a larger sub- volume and 
at large N is needed [45j]. Owing to the oscillatory be- 
haviour of the integrand, the methods presented here, 
numerical integration or the saddle point approximation, 
do not provide access to this parameter range. The os- 
cillatory integrand causes problems with the numerical 
integration of Eq. (|15[) while the rather complicated an- 
alytic structure of the grand canonical partition function 
at complex /i interferes with the saddle point approxima- 
tion. 



IV. CONCLUSIONS 

We have explored the properties of the probability dis- 
tributions P{N) of the net baryon number in systems 
that exhibit a second order chiral phase transition. The 
critical behavior was modelled within the Landau theory 
of phase transitions, where the effective thermodynamic 
potential is expressed as a polynomial in the order pa- 
rameter. The coefficients of the effective potential were 



parametrized such that the scaling behavior of the spe- 
cific heat and of the cumulants of the net baryon number 
in the mean-field approach and in the 0(4) universality 
class, respectively, is reproduced. The resulting effective 
potential preserves the baryon-anti baryon symmetry and 
the periodicity in the imaginary baryon-chemical poten- 
tial, obtained in QCD. 

In the mean-field approach, where the specific heat 
critical exponent a — 0, the probability distribution of 
the net baryon number was computed analytically. In 
the model which reproduces the 0(4) scaling of the net 
charge fluctuations, the probability distribution could 
not be obtained in closed form. In this case we applied 
the numerical and the saddle point methods to calculate 
P(N) and we discussed the validity and the limitations 
of these methods. 

We have shown that in the mean-field case, the contri- 
bution of the singular part of the thermodynamic poten- 
tial, which is responsible for the critical behaviour of the 
order parameter, results in a broadening of the probabil- 
ity distribution. On the other hand, in the model, which 
exhibits 0(4) scaling of the net baryon number fluctua- 
tions, the contribution of the singular part results in the 
opposite behaviour, i.e. in a narrower net baryon number 
probability distribution. However, using the numerically 
computed 0(4) probability distribution P(N), we are not 
able to reproduce the critical behaviour of the net baryon 
number cumulants. The main reason for this failure is, 
that the resulting probability distribution P(N) is appli- 
cable only up to a finite \N\ < 50 and in a rather small 
volume. In order to reproduce the singular structure of 
higher-order cumulants directly from the probability dis- 
tribution one would need to compute P(N) for large net 
baryon number N and for a sufficiently large volume. 
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